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DEPOLARIZATION CANALS AND INTERSTELLAR 
TURBULENCE 

Andrew Fletcher 1 and Anvar Shukurov 1 



Abstract. Recent radio polarization observations have revealed a ple- 
thora of unexpected features in the polarized Galactic radio background 
that arise from propagation effects in the random (turbulent) interstel- 
lar medium. The canals are especially striking among them, a random 
network of very dark, narrow regions clearly visible in many directions 
against a bright polarized Galactic synchrotron background. There are 
no obvious physical structures in the ISM that may have caused the 
canals, and so they have been called Faraday ghosts. They evidently 
carry information about interstellar turbulence but only now is it be- 
coming clear how this information can be extracted. Two theories for 
the origin of the canals have been proposed; both attribute the canals to 
Faraday rotation, but one invokes strong gradients in Faraday rotation 
in the sky plane (specifically, in a foreground Faraday screen) and the 
other only relies on line-of-sight effects (differential Faraday rotation). 
In this review we discuss the physical nature of the canals and how 
they can be used to explore statistical properties of interstellar turbu- 
lence. This opens studies of magnetized interstellar turbulence to new 
methods of analysis, such as contour statistics and related techniques 
of computational geometry and topology. In particular, we can hope 
to measure such elusive quantities as the Taylor microscale and the 
effective magnetic Reynolds number of interstellar MHD turbulence. 



1 The promise and pitfalls of radio polarization 

High-resolution maps of the Milky Way synchrotron emission show a complex, 
irregular system of polarized and depolarized structures (Wieringa et al. 119931 
Uyamker et q/. ll998bl Duncan et al- lRKM Gray et a/. [1999, Haverkorn et aZ. IMTUl 
Gaensler et al. 121)1)11 Wolleben et al. I2l)l)6f) . Many polarized structures can be 
attributed to objects in the interstellar medium (ISM), such as supernova rem- 
nants and shock fronts. These objects are usually bright in total intensity as well. 
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Fig. 1: A 24° x 9° section of the 1.4 GHz Effelsberg Medium Latitude Survey cen- 
tred at (l,b) = (162°, 0°) (Reich et al. ^M^ . Top: total intensity, with large-scale 
structure added from Dwingeloo data. Middle: polarized intensity, also including 
the Dwingeloo large-scale structure. Bottom: polarized intensity observed with 
the Effelsberg telescope with large scale structure missing. Note the strong dif- 
ference between the (preliminary) absolutely calibrated polarized intensity map 
[middle) and the map with artificial base-levels {bottom). 



However extended, diffuse polarized emission is also rich in structure which most 
plausibly arises from propagation effects as the radiation passes through the ran- 
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dom ISM. As illustrated in Fig. ^ much of this structure is not mirrored in the 
total synchrotron intensity so the structure in polarization is not primarily caused 
by variations in emissivity. At low frequencies, the cloudy appearance of polar- 
ized intensity maps differs between adjacent, narrow frequency bands (Haverkorn 
et al. 2003) strongly suggesting that the polarization varies with wavelength and 
so a connection with Faraday rotation is plausible. Random fluctuations in the 
magnetized ISM — caused by its turbulent flows, for example — leave their fin- 
gerprint in the random appearance of polarized intensity maps. If a connection 
can be established between measurable properties of the radio maps to the physics 
of the ISM, a powerful new source of information about the dynamic state of the 
interstellar plasma will become available, particularly its little studied magnetic 
properties. Similarly, a reliable separation of the Milky Way synchrotron polar- 
ized foreground from cosmological signals, either through the creation of templates 
from low-frequency data or through statistical analysis, requires a thorough un- 
derstanding of the origins of structure in the diffuse polarized emission. In the 
context of cosmological studies, the effects of the ISM on the polarized emission 
propagating through it are especially important. 

One of the most striking features of the maps is the twisting network of narrow, 
dark canals running through regions of bright polarized intensity, clearly visible 
in Fig. ^ at 1.4 GHz and common in other maps at this frequency (Uyamker et al. 
199813 Gaensler et al. [2001) and longer wavelengths (Haverkorn et al. 2000). The 
observed canals have the following properties: 

1. the observed polarized intensity falls to zero, because the polarization angle 
changes by 90° across the canal; 

2. a canal is one telescope beam wide; 

3. a canal passes through a region of significant polarized intensity; 

4. the canal is not related to any obvious structure in the total intensity. 

Care in interpreting the observations is especially important in the case of 
polarized emission: interferometer observations miss emission at large scales (de- 
termined by their shortest baseline) , and single-dish surveys often have arbitrarily 
set base-levels around their edges. Correcting for the missing large-scale structure 
involves adding a smoothly varying signal to the Stokes parameters. In the case 
of total intensity (Stokes /) the consequences are easily predictable. However, 
the Stokes parameters Q and U are not positive definite and therefore algebraic 
addition of a smooth component can create and/or remove zeros in either. Hence 
minima in polarized intensity (Q 2 + C/ 2 ) 1 / 2 can become maxima and vice versa. 
The effect of absolute calibration on the polarization angle \P = | arctan (U/Q) 
can be equally counter-intuitive. Reich et al. I|200 4} clearly identify and discuss 
these problems. The dramatic effect the missing data can have on the appearance 
of the polarization pattern can be seen by comparing the bottom two panels of 
Fig. G] It is clear that unless statistical parameters can be identified that are 
independent of the Q and U intensity values, polarization surveys should be tied 
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to an absolutely calibrated reference frame (Uyaniker et al. 19983 Reich et al. 
2004). The recent absolutely calibrated all-sky survey at 1.4 GHz by Wolleben et 
al. (J5006 ) will be a mine of useful information for polarization studies. Haverkorn 
et al. ( 2005 argue that if fluctuations in Faraday rotation measure are sufficiently 
strong and of a small scale, the resulting structure in Q and U can be of such 
a small scale that the missing short spacings will not cause a problem. In other 
words, Q and U may have no large scale structure; this can easily happen at metre 
wavelengths. 

Two theories accounting for the origin of the canals have been proposed: one in- 
vokes strong plane-of-sky gradients or discontinuities in Faraday rotation measure 
in a foreground screen between the emitting region and the observer (Haverkorn 
et al. 200U1 12U04I Haverkorn & Heitsch 2004 ) and the other proposes differential 
Faraday rotation along the line-of-sight through the emitting layer (Beck 1999 
Shukurov & BerkhuijsenEIIIi- 

In this review we summarize recent work on the interpretation of depolarized 
canals and provide a more detailed explanation of some analytic results than is 
available in the readily accessible literature. After introducing basic equations and 
defining our notation in Section |3 we describe the effects of differential Faraday 
rotation and Faraday screens, demonstrate how canals are formed by the two pro- 
posed mechanisms and discuss how observable quantities behave in the vicinity of 
a canal in each case. Then in Section|3we turn to the interpretation of the proper- 
ties of the canals in terms of the physical state of the ISM. In Section ITTI the case 
of differential Faraday rotation is discussed in the context of the statistical prop- 
erties of contours of a random field and we show how canals of this type can yield 
information about the smallest scales of ISM turbulence. Similar techniques can 
be applied to canals produced by Faraday screens, and more generally to contours 
of any randomly distributed observable quantity. Section 12.21 demonstrates that 
true discontinuities in Faraday rotation are required if a foreground screen is to 
produce canals similar to those observed. In the case that these discontinuities are 
shock fronts originating from supernovae, we show in Section [3. 21 how the distance 
between canals is related to the separation of shock fronts. 



2 The complex polarization, Stokes parameters and depolarization 

The complex polarization is commonly written as 

/ v W/(r)e(r)exp[2»^(r)]d^ 
V - P ° f v W(r)e(r) dV ' (2J) 

where po — 0.7 is the maximum degree of polarization for synchrotron emission 
(a function of the spectral index), W is the beam profile around a given position 
in the sky, e is the synchrotron emissivity, and ip is the local polarization angle. 
The integrals are taken over the volume of the beam cylinder and r = (x, y, z) 
specifies a location with respect to the sky plane (x,y) and the line of sight, z. 
The modulus and argument of V are the observed degree of polarization p and 
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polarization angle \E' respectively, 

V = p exp (2i\£). (2.2) 

When polarized radio emission propagates through magnetized and ionized ISM, 
the local polarization angle -0 undergoes Faraday rotation: 

0(r) = o (r) + 0(r) , (2.3) 

where tpo is the intrinsic polarization angle (perpendicular to the magnetic field 
component transverse to the line of sight at the point of emission) , 

/>OC 

<t> = KX 2 n e (z')B z (z')dz', (2.4) 



K = 0.81 radm~ 2 pc _1 cm 3 /xG _1 is a constant, n c is the density of free thermal 
electrons, B z is the component of the magnetic field along the line of sight z, A is 
the wavelength and the observer sits at z — oo. The Faraday depth, at a position 
(x, y) in the plane of the sky, is defined as, 



F{x,y) = <j>(-oo) =K\* / n c {z)B z {z)dz . 

and <fi(z) is known as the Faraday depth to a position z (cf. Eilek 1989). F and 
4>{z) are convenient variables to work with; <p(z) gives the change in polarization 
angle of a photon of wavelength A passing through the region z < z' < oo, and 
F is the maximum amount of Faraday rotation in a given direction in the sky. 
Note that F/X 2 can differ strongly from the observed Faraday rotation measure 
RM = difi/dX 2 . For example, RM = F/2X 2 when emission and rotation occur 
together in a uniform slab (see below). However, RM = F/X 2 for a layer that only 
rotates — a Faraday screen. A more comprehensive discussion of depolarization 
and the effects of Faraday rotation can be found in, e.g., Burn l|1966f) and Sokoloff 
et al. lfT§9"g|) . 

Observations of linearly polarized emission measure the Stokes parameters J, 
Q, U and these are related to V, p and "J: 



v = ^±, p = ^ , *=i 



Arctan — — |7r(signQ — l)signf/ 



iu VQ 2 + u 2 
—> p = T 

(2.5) 

where | Arctan x \ < ir/2 denotes the principal value; the above form of ^ allows 
for a consistent choice of the branch of arctan for any combination of signs of Q 
and U, so that -tt/2 < * < vr/2. 

2.1 Differential Faraday rotation 

If linearly polarized emission undergoes Faraday rotation, the photons originating 
at different positions along a single line of sight will suffer different rotations of 
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Fig. 2: Diagram showing how canals are formed by differential Faraday rotation. 
The observer is on the right and x is oriented perpendicular to their line of sight; 
each horizontal strip represents one telescope beam. The Faraday depth varies in 
the plane of the sky (with x) as shown with thick solid line; the intrinsic polar- 
ization angle is uniform (vertical double arrows on the left). Emission generated 
on the x-axis cancels the emission generated in a layer where the former is ro- 
tated by 7r/2, i.e., where Fix) = tt/2. Layers whose emissions mutually cancel are 
indicated by perpendicularly hatched boxes. Only the emission from unshaded 
layers reaches the observer. The arrows on the right show the observed polar- 
ization plane which has been Faraday rotated in layers to the right of the visible 
(unshaded) layer to give the polarization angle ^> in the range indicated. Note 
that layers to the right of the 'source' (an unshaded layer) act as a Faraday screen, 
and so the rotation angle is equal to F. The canals occur where the polarized 
emission is fully cancelled (regions labelled B). The emission observed at values 
of F(x) slightly smaller than nir (regions A) and slightly larger than nir (regions 
C) originates in different layers. The difference in F(x) between regions A and 
C is 7r/2. Therefore, emission from region A is rotated by tt/2 before reaching 
the observer, whereas emission from region C is rotated only slightly. Thus, the 
polarization angles observed on the two sides of a canal differ by tt/2. A canal 
that occurs when F(x) = nir with n = 1 is shown near the top, and with n = 2, 
near the bottom. (Adapted from Shukurov & Berkhuijsen 2003 ) 
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(a) Differential Faraday rota- 
tion 



(b) Faraday screen 



Fig. 3: The relative positions, along the line of sight, of the regions producing 
synchrotron emission and Faraday rotation. The observer is on the right, (a): the 
case of differential Faraday rotation, i.e., both emission and rotation occur in the 
same region, (b): a Faraday screen, where all the emission is produced behind 
the rotating layer. 



their polarization planes. So even if all the photons are emitted with the same 
polarization plane, and the amount of Faraday rotation per unit path length is 
constant, the emerging radiation will comprise a mixture of polarization angles 
and the resulting degree of polarization will be reduced. This process is called dif- 
ferential Faraday rotation (Fig. |2J) ; the depolarizing effect is wavelength-dependent 
but occurs even for infinitely small telescope beams. 

Although the effect appears rather simple to understand, one can discover un- 
expected subtleties from simple mathematical analysis. Consider a slab extending 
between —h<z<h along the line of sight, with constant synchrotron emissiv- 
ity, e = const, and uniform magneto-ionic medium, n c B z = const, throughout as 
shown in Fig. |3(a)| The telescope beam is assumed to be infinitely narrow for 
simplicity, i.e. W is a ^-function in x, y; this makes trivial integration in the sky 
plane in Eq. (|2.ip . Equation l|2.4|) then yields 



where F = 2KX 2 n e B z is the Faraday depth of the slab. Then, for ipo — const, 
Eqs. and (|2.3[> lead to the well known result (Burn ll9bl)J) 




z < —h , 
\z\<h, 
z > h , 



(2.6) 



sin F 



exp[2i (V>o + i-F)]. 



(2.7) 



F 



To see more clearly how differential Faraday rotation produces depolarization 
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(a) Differential Faraday rotation 



(b) Faraday screen 



Fig. 4: Variation of Q (solid) and U (dashed) (in arbitrary units) across depo- 
larization canals caused by differential Faraday rotation (a) and in a nonuniform 
Faraday screen (b). In each panel, the canal occurs at x — X c = 0.5 where 
Q = U = 0. In Panel (a), the beam is infinitely narrow and the Faraday depth 
varies linearly with x smoothly passing through F — ir at x — 0.5. For the Fara- 
day screen, the beam is a Gaussian of FWHM D = 0.2, with the Faraday depth 
changing discontinuously by AF = at x — 0.5. The Stokes parameters Q and 
U were obtained from Eqs. H2.1[) - I|2.5[l for the respective configurations shown in 
Fig. |3 with e, n c B z and ipo uniform along the line of sight. 



canals, it is useful to rewrite this with the help of Eqs. (|2.2I) and Q2.5[l : 



9. 
i 

u 
7 



Po 



Pa 



sinF 



F 
sin_F 



cos2(V'o + 5 F )' 



F 

Po— —sm2(ip + lF). 



(2.8) 
(2.9) 
(2.10) 



Equation H2.8(l shows that p = where \F\ = tvk, n = 1, 2 If F is a continuous 

function of position with well-behaved contours, canals will develop on the contours 
of F where \F\ = nn. Therefore, for a canal produced by differential rotation we 
know the magnitude of the Faraday depth along the canal (up to an integer factor 
n) . Techniques for retrieving further information about the physical state of the 
ISM from the statistical properties of these contours is discussed in Section 

We are free to measure the polarization angle from a line perpendicular to the 
transvese magnetic field near the canal; then ip — 0. Equations (|2.9I) and (|2.10|l 



A. Fletcher & A. Shukurov: Depolarization canals 



show that in this reference frame 

Q oc sin F cos F , U oc sin 2 F . 

When F{x) passes through an integer multiple of 7r, both Q and U vanish to 
produce a canal, but Q will change sign across the canal and U will not. In other 
words, there is always a reference frame in the sky plane such that one of the 
Stokes parameters changes sign, but the other does not across a canal produced 
by differential Faraday rotation. This behaviour is illustrated in Fig. 4(a) and we 



will see in Section that this does not occur in the case of canals formed in a 



Faraday screen as illustrated in Fig. 4(b) 



For the signs of Q and U shown in Fig. 4(a) W is in the second quadrant to the 
left of the canal and in the first quadrant on the right of it. As Q — > and U — > 
for x — > X c , where X c is the position of the canal, we have W = — 7r/2 for a; = X c — 
but W = 7r/2 for x = X c + — see also Eq. 1)2. 5[l . Thus the polarization angle 
changes discontinuously across the canal although the magneto-ionic medium is 
perfectly continuous. In particular, RM has the same values on both sides of a 
canal (but is undefined at it). As illustrated in Fig. [3 the jump of the polarization 
angle in a continuous medium occurs because the polarized photons detected on 
the either side of the canal originate at different depths in the source and therefore 
have experienced different amounts of Faraday rotation. 

The discussion above refers to the case of an infinitely narrow telescope beam; 
these canals are infinitely narrow curves in the sky plane. The effect of a finite 
beam width is twofold. Firstly, the canals in observed maps will be one beam 
wide, as the telescope combines emission from either side of the contour where 
\F\ = wk. Secondly the canals will gradually fill up with polarized emission as 
the resolution coarsens, so that asymmetry of the magneto-ionic medium within 
the beam becomes noticeable. 

Canals are also filled up if the synchrotron cmissivity varies along the line of 
sight provided this variation is not compensated by a variation in Faraday rotation. 
It can be shown, in particular, that a double-layered source where synchrotron 
intensities produced in each layer are I± and I2, with \I± —I^l/Ii = s, and F = nit, 
will produce a partially filled canal with the fractional polarization p/po — e/n, 
valid for e <C 1, i.e., for weak inhomogeneity. This effect can explain why the 
observed canals are not all closed curves as might be expected of contours of a 
continuous function. 

Tangled magnetic fields and fluctuations of the electron number density can 
produce additional depolarization known as internal Faraday dispersion, charac- 
terized by the standard deviation of the Faraday depth, of- This effect also fills up 
a canal, so that the minimum fractional polarization is p/po — o F /n for op <C 7r. 
The above two effects are discussed by Fletcher & Shukurov (in preparation). 

Joint action of various depolarization effects can produce counter-intuitive ef- 
fects. For example, Sokoloff et al. (1998- their Sect. 6.2 and Fig. 5 in Erratum) 
show that internal Faraday dispersion combined with differential Faraday rotation 
can lead to a jump of 90° in polarization angle whilst the degree of polarization 
varies smoothly and p differs from zero (if observed with infinitely narrow beam) ; 
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Fig. 5: Sketch showing the variation of Faraday depth F across the sky, here along 
x. There is a constant gradient of F in the range xo < x < xo + Ax, and the 
telescope beam has a width D. 



this would appear as a canal when observed at finite resolution because of the 
jump in the polarization angle. 

The position of a canal produced by differential Faraday rotation must change 
with the wavelength at which it is observed because F depends on A according to 
Eq. (|2.4(l . To a first approximation, a change of wavelength from A to A' = A + SX 
will displace the canal by Sx where 



Sx 



SF 



NF 



SF = F(A') - F(A), |<5A|<A 



For example, if differential Faraday rotation is the only depolarization mechanism, 
we have F = 2RM A 2 with RM being independent of A, and then 



Sx ~ 2L 



RM 



5X 

T 



RM 



VRM 



where Lrm is the scale at which RM varies. Shukurov & Bcrkhuijsen (2003) 
argue, using a similar estimate, that the displacement can be difficult to observe 
if ^ -V which is often the case. 



2.2 Foreground Faraday screens 

If synchrotron emission and Faraday rotation occur in distinct layers — with the 
rotating layer nearer to the observer as illustrated in Fig. 3(b) — the foreground 



magneto-ionic medium is called a Faraday screen. If Faraday rotation in the screen 
is uniform across a telescope beam then no depolarization occurs; the polarization 
angle of radiation emerging from the emitting layer is rotated by A^> = F, but 
the mixing of different angles within a beam, that causes depolarization, does not 
occur. A Faraday screen can only cause depolarization if it is nonuniform on the 
scale of the beam, so that different lines of sight within the beam undergo differing 
amounts of Faraday rotation. 
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Fig. 6: The variation of the fractional polarization and polarization angle across 
a Faraday screen with a smooth variation of Faraday depth described in Fig. 
with xq — 0.4, Aa; = 0.2 and X c — 0.5. Left-hand panels: (a) the change in 
the beam-averaged degree of polarization (solid) and polarization angle (dashed) 
for AF = ir/2 convolved with a square beam of width D = 0.2. The fractional 
polarization does not reduce to zero although $ changes by ir/2 across the region. 
The reason is illustrated in panel (b) where we show the change in the polarization 
angles obtained along individual lines of sight (dashes) within the beam (circle). 
It is clear that the polarization cannot drop to zero in this configuration. The 
right-hand panels show the same as on the left, except that here AF = it. Now, 
panel (c) shows that the fractional polarization does drop to zero at x = X c 
producing a canal, but the corresponding change in 'J' across the canal is n. Panel 
(d) explains how this variation in ^ among individual lines of sight leads to the 
complete depolarization. The calculation used Eqs. I2.1H2.5I for the configuration 
of Fig. El 



Consider a uniform polarized radio source shining through a Faraday rotating 
layer. Consider a simple variation of the Faraday depth in the sky plane shown 
in Fig. where F changes by AF across a distance Ax, and the width of the 
telescope beam is D. We assume that the telescope beam has a square profile for 
simplicity. In the regions where F = const in Fig. Eq. (|2.1|l reduces to 

f exp[2i(V>o + F )] , x<x -D/2, 

\ exp[2i(V>o + F + AF)] , x > x + Ax + D/2 , [ ' 

so p — KcV = po and no depolarization occurs if F is constant. In the region 
xo + D/2 < x < xo — D/2 + Ax, where F varies linearly, Eq. (|2.1|l yields 

1 / x -xn\ 

V = po exp[2i(V> + F )]- / exp 2iAF D — dx 

U Jx-D/2 V U J 



12 



Title : will be set by the publisher 



= Po^^exp[2^o + F)], (2.12) 

where AFd = D dF/ dx is the increment in the Faraday depth across the telescope 
beam arising from a continuous gradient in F. (We ignore, for simplicity, the 
ranges \x — xq\ < D/2 and \x — xo + Ax\ < D/2 with a mixture of constant F and 
a gradient in F within a beam.) Then 

sin AF D 



P = Po 



AF 



D 



(2.13) 



Thus, the depolarization in a continuously inhomogeneous screen is complete if 
|AFd| = mr and we can then expect a canal to form along a contour where 
|VF| = nir/D (where we have replaced dF/dx by VF). An increment of |AFd| = 
mr across the beam means that the change in polarization angle across the canal 
is A^ = 180°, i.e., no change — see Ea. 12.121 However, canals attributed to a 
Faraday screen by Haverkorn et al. (2000) have AW = 90°. As we illustrate in 
Fig. this change in W in fact cannot produce a canal if F varies continuously: 
a gradient in F that produces a 90° change in polarization angle across the beam 
does not lead to complete depolarization, whereas a gradient that produces 180° 
of Faraday rotation does give p = 0. Canals with AW = (or 180°) predicted here 
have not yet been observed. 

Depolarization can be complete if the change in Faraday rotation within the 
beam produces a mixture of polarization angles that exactly cancel each other. 
Complete cancellation depends on both the variations in the Faraday depth across 
the screen AF and the telescope resolution, so that the pattern of canals can 
change very significantly if observed at different resolution. 

Figure [7| shows how the depth of a canal produced by a Faraday screen depends 
on the increment in Faraday depth across the beam AFp and on how well resolved 
this increment is, Ax/ D. We can see that all gradients in F can produce p = if 
observed at an appropriate resolution and so we should expect to see canals with 
all jumps in polarization angle across them. The only condition for that is that the 
total variation in F available exceeds 7r/2, the minimum that can produce a canal 
for Ax/D <C 1. The jump in W across the canal (which is one beam wide) in Fig.0 
is equal to the corresponding value of AF (shown in the legend) if Ax/D < 1 and 
to AFd otherwise. 

It is evident from Fig. that the only way that canals with AW = 90° will be 
preferentially created is if the variations in F in the Faraday screen are discontin- 
uous. In this case complete cancellation will only occur if all of the polarization 
angles within a beam are a mixture of exactly orthogonal pairs and canals will be 
formed along contours where AFd = (2n + \)n/2. Like the case of differential 
Faraday rotation, canals produced by discontinuities in a Faraday screen do not 
depend on the resolution of the observations. In Section 13.21 we will discuss how 
the statistical properties of canals formed by discontinuities in a Faraday screen 
can provide information about the shock-wave turbulence in the diffuse ISM. 

For a canal produced by a discontinuous change of AF, the behaviour of Q 
and U across the canal can obtained from Eq. H2.11J1 with AF = (2n + 1)tt/2: each 
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Fig. 7: The minimum degree of polarization produced by a gradient of the Faraday 
depth in a foreground screen, p m i n , as a function of the relative extent of the 
gradient Ax/D, where D is the width of a square beam. Different curves represent 
different Faraday rotation increments (AF, shown in legend). For Ax/D < 1, F 
varies over a range narrower than the beam width, and is effectively discontinuous 
when Ax/D <C 1. Complete depolarization in a discontinuous screen, Ax/D — > 0, 
occurs if AF = (2n + 1)tt/2, n = 0, 1, . . .. 



changes by n (neglecting 2mr in the arguments of cos and sin), 



and likewise for U. Figure 01 shows an example of how Q and U change across 
a canal produced in a Faraday screen and by differential Faraday rotation. In 
contrast to the case of differential Faraday rotation, Eqs. I|2.9[) and (|2.10() . both Q 
and U of Eq. 12.141 will always change sign across a canal produced by a Faraday 
screen. (A special case is where Q (or U) is zero whereas U (or Q) changes sign at 
the canal.) This offers an opportunity to reveal the nature of an observed canal 
from a careful analysis of the variation of Q and U across it. 

3 What can we learn about the ISM from depolarization canals? 

If the mechanism which forms a canal can be identified — for example by looking 
at the variation of Q and U across the canal — we then know the location of level 
lines of Faraday depth where \F\ — nir, in the case of differential rotation, or the 
projected position of discontinuities with |AFd| = (2n+ l)%/2, for the action of a 
Faraday screen. In both cases, canals represent contours of a random function of 
position in the plane of the sky, namely the Faraday depth or the magnitude of its 




x < X c 
x > X c 



(2.14) 
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Fig. 8: Crossings of a level F c by a random function F(x) are indicated with 
vertical dashed lines. 



gradient. In this section we discuss how statistical properties of the contours of a 
random function can be used to measure its parameters, and hence how statistical 
properties of the depolarization canals can be interpreted in terms of the properties 
of the random (turbulent) ISM. 

3.1 Differential Faraday rotation: contour statistics and turbulence 

As we show here, the mean separation of the contours of a differentiable random 
function is sensitive to the curvature of its autocorrelation function at zero lag, 
i.e., to a parameter known as the Taylor microscale in turbulence theory. Reviews 
of this theory can be found in Longuet-Higgins l|19fi7J) . Sveshnikov (1966) and 
Vanmarcke (1983). Observations of canals (or contours of another observable), 
obtained under quite a modest resolution, provide a possibility to measure the 
Taylor microscale, related to the dissipation scale of the turbulent flow and to the 
Reynolds number. Contour statistics have been used in astrophysics in the studies 
of structure formation (Peebles 119841 Barden et dl. 1986 Ryden 1988; Ryden et 
al. 1989) and the cosmic microwave background (Coles & Barrow 1987). 

3.1.1 Overshoots of a random function in one dimension 

We start with crossings of a random function of a single variable over a refer- 
ence level (Fig. IHJ); a clear and detailed discussion of this problem can be found 
in §9 of Sveshnikov ( 1 966fl . Consider a stationary, differentiable random function 
of position J 7 (x), whose realisation, or measured value, is denoted by F(x). (We 
have in mind applications where F can be the Faraday depth or the magnitude 
of its gradient.) The function will pass, from below, through the level F c in the 
infinitesimal interval (x, x + dx) if the following inequalities are satisfied simulta- 
neously: F(x) < F c and T(x + dx) > F c . For a continuous ^(x), we can write 
T(x + dx) w T(x) + J-'(x) da;, where T' = dT /dx. Thus we have 



F c - T'{x) dx < T{x) < F c for T'(x) > 



(3.1) 
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In order to calculate the probability that these two inequalities hold, consider the 
joint probability density of T and its derivative T' at a given position x, which we 
denote f(F,F'). Then the probability of a overshoot above a level F c is given by 



poo nF c 

P{F C - T'(x) dx < T{x) <F C }= / AF" / dF f(F, F' 

JQ JF c -F'(x) dx 



where integrations are performed over all values of F and T' that satisfy Eq. (|3.1|1 . 
The integral over F is immediately taken using the mean value theorem: 



/ " f(F,F')dx = f(F c ,F')F'dx, 

J F c -F'(x)dF 



and then 



/•OO 

P{F C - T'{x) dx < T{x) < F c } = dx f(F c , F')F' dF' . (3.2) 

Jo 

It is convenient to introduce the probability density of an overshoot, i.e., the 
probability p a (F c ) of an overshoot above a level F c at a position x per unit interval 
of x. The above result can be written as 



/ f(F c ,F')F'dF' . 
Jo 



Similarly, we obtain the probability density for F(x) to cross the level F c from 
above: 



Pb(Fc) = - I f(F c ,F')F'dF' , 

J — oo 



and the total probability density for F(x) to cross the level F c at a position x 
either from above or from below is given by 



Pc 



/CO 
f(F c ,F')\F'\dF' . (3.3) 
-oo 



It can be shown (Sveshnikov I1966f) that the mean numbers of crossings of a 
level F c per unit length, A^, is equal to p c , i.e., the mean separation of the crossing 
points is given by 

£c = N\ X =P~ X , 

where overbar denotes averaging. 

The results can be written in a closed, explicit form if T[x) is a Gaussian 
random function because then the fact that T(x) and F'(x) are uncorrelated 
(which is true for any stationary, differentiable random function) implies that 
they are statistically independent, so that 

f(F,F') = f F (F)f F ,(F'), 
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and hence (we note that F' = 0) 
1 



f(F,F') 



a F (2iry/ 2 



exp 



2a\ 



a F >(2ny/ 2 



exp 



(F' 



2a F , 



where we have introduced the standard deviations of F and i 7 ", a F and cr F /, 
respectively. Integration over F' in Eq. (|3.3|) can now be easily performed to 
obtain 



AT 1 ° F 

Ni = p c = exp 

7r a F i 



(F C -FY 



2a F 



(3.4) 



In application to canals, where F c = mr, n = 1,2, ... in the case of differential 
Faraday rotation, one needs to calculate the mean separation of positions where 
|-F(x)| = F c . The corresponding probability density involves an additional contri- 
bution arising from crossings at F(x) — —F C) which should be added to the right- 
hand side of Eq. 1)3.4(1 . The additional term is proportional to exp[— (F c +F) 2 /2a F }. 
This term can be significant if F « 0, so that \F C — F\ ~ \F C + F\. If only one of 
the two terms is significant, it is sufficient to replace F by \F\ in Eq. (|3.4|l . 
Since 



d 2 C(r) 



dr 2 



where C(r) is the autocorrelation function of F(x) normalised to C(0) = 1 and 
It is known as the Taylor microscale (e.g., Sect. 6.4 in Tennekes & Lumlev fl972f) . 
we obtain the following expression for the mean separation of the positions where 
\F(x)\=F c : 



= —= It exp 
v2 



2a F 



(3.5) 



3.1.2 Contour statistics in two dimensions 

The theory of overshoots can be generalized to two dimensions, where we consider 
a random field !F(x) depending on two variables x = (x,y) (e.g., position in the 
sky plane); we also introduce V = VJ 7 . The probability of an overshoot is now 
equal to the probability that the following inequalities are satisfied simultaneously: 

^(x) < F c , .F(x') > F c , x' = x + dx , 

where dx = Vdx, with V = VJ-/\VT\ the unit vector parallel to the gradient of 
T . Expanding .F(x') ~ ^(x) + VJ 7 ■ dx and retaining only the term linear in dx, 
we reduce the above inequalities to 

F c -|V^(x)|da;<^(x)<^ c , 

and, using arguments similar to those that lead to Eq. I|3.2(l and further, we obtain 
the probability of a crossing (either upward or downward) of the contour ^F(x) = F c 
in two dimensions: 



/oo />oo 
/ /(F c ,V)|V|d^dn 
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where f(F c , V) is the joint probability distribution of T and its gradient. This 
quantity can be identified with the mean length of the F c contour per unit area, 

N 2 = P W . 

In both one and two dimensions, the probability of crossing a contour is given 
by the mean value of the magnitude of the gradient of the random field at that 
contour. As discussed by Ryden l|1988[) . for an isotropic random field F, the mean 
number of crossings per unit length along a line, Ni, and the mean length of a 
contour per unit area, N2, are related, 

N 2 = , 

and so statistical properties of contours in the plane can be conveniently studied 
using one-dimensional cuts, which is technically easier. 

3.1.3 The separation of canals and parameters of ISM turbulence 

Shukurov & Berkhuijsen (2003) applied a one-dimensional version of this theory 
to canals observed at A20.5cm in the direction to the Andromeda nebula (around 
I = 120°, b = —20°). Their estimate of the mean separation of canals is 5' which 
results in a estimate of the Taylor microscale of RM (or, more precisely, the Faraday 
depth) of It ~ 2'. For the presumed mean distance within the Milky Way's non- 
thermal disc in that direction, 1 kpc, the corresponding linear scale is 0.6 pc. If this 
scale can be identified with the Taylor microscale of interstellar turbulence (which 
may or may not be the case!), the effective value of the Reynolds number follows 
as Re ~ (Io/It) 2 — 10 5 for the energy-range scale of Zo = 100 pc. This value of Re 
is many orders of magnitude smaller than those based on the molecular viscosity. 
If this estimate can be confirmed by more comprehensive and deeper analyses, 
this would imply that the dissipation range of the interstellar MHD turbulence is 
controlled by relatively large viscosities, perhaps associated with kinetic plasma 
effects. The latter are often described as anomalous viscosity or resistivity; useful 
recent results were obtained in the context of galaxy clusters by Schekochihin et 

al. CJE3- 

3.2 Faraday screens: the mean separation of canals and shocks 

If canals are caused by rotation measure discontinuities in a Faraday screen, the 
mean distance between them will be related to the distribution of discontinuities 
in the ISM. Most plausibly these are shock fronts. (Another possible source are 
tangential discontinuities but it is difficult to identify a physical reason for their 
widespread occurrence.) 

A canal forms where the background Faraday depth F changes by |AF| = 
(n + 1/2)tt across a shock front. If the magnetic field is frozen into the gas and 
the gas density increases by a factor e p then so will the field strength (we neglect 
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a factor ~ 1/ \/2 due to the relative alignment of the shock and field) . Therefore 

L. (3.6) 



AF 



The increase in gas density depends on the Mach number of the shock M. as (e.g. 
Landau & Lifshitz ITMIjl 



{j + l)M 2 
(7- l)M 2 + 2 



M 2 



M 2 



(3.7) 



where we have used 7 ~ 5/3 for the polytropic index of the gas. Combining 
Eqs. I|3.6[) and (|3.7|l . we obtain the minimum Mach number necessary to produce 
a canal, 




3(\AF/Fo\ + iy/ 2 



-1/2 



7T 

2' 



(3.8) 



As follows from Sect. 12.21 stronger shocks can also produce canals. 

The dominant source of the primary shocks are supernovae. When these shocks 
encounter gas clouds they will be reflected, generating a population of secondary 
shocks. Following Bykov & Toptygin l)1987|l . the probability density function of 
both primary and secondary shocks is given by 

1 , 3C(a)/ c i 



G(M) ~ Go 



(M - I) 4 



(3.9) 



where Go = S^tttq, S is the supernova rate per unit volume, r$ is the maximum 
radius of a supernova (primary) shock, and a characterizes the dependence of the 
shock radius on the Mach number, r{M) — roAi~ a / 3 , f c \ is the volume filling 
factor of diffuse clouds and C(a) is a constant. For the Sedov solution, a = 2 and 
for a three-phase ISM, a = 4.5; C(a) = 2.3 x 1CT 2 for a = 2 and 4.1 x 1(T 3 for 
a = 4.5. The mean number of shocks of a strength > M. crossing a position in 
the ISM per unit time is 



H(M) 



G{M)dM = Go 



M 



1 



C(a)f c 



aM c 



(M - l) 3 



and the mean separation of shocks with a strength > M. is 

L(M) 



(3.10) 



(3.11) 



H{MY 

with c s the sound speed. Substituting Eqs. I|3.8|l and l|3.10|l into Eq. l|3.11|l . we 



find the mean separation of shocks in three dimensions: 

. -1 / r, \ 2 

L ~ 9pcf' 



X 



.lOkms 
1 



0.02 yr 



,•-1 



R 



15kpcy V50pc7 ViOOpc 



ro 



C(a)f, 



cl 



olM* (M, - l) 3 



(3.12) 
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where ^sn is the Galactic supernova rate, R and h are the radius and scale height 
of the star-forming disc respectively; the term in the square brackets is close 0.2 
for M. — 1.2 and a = 4.5. Projecting onto the sky plane changes £ by a factor 
£ = 7T _1 J w sin 2 9d9 = 1/2. If the shocks fill a region of a depth d, the separation 
of the shocks in the plane of the sky is further reduced by a factor d/L. The 
angular separation C at an average distance d/2 to the shocks follows as 

C = * (3.13) 

As an example, let us use the observations by Haverkorn et al. ( 2003) as inputs 
to Eq. I|3.13l) and compare the observed separation of canals with the predicted 
separation of shocks. In this held the average Faraday depth is F ~ tt radians and 
the most abundant canals are produced when \AF\ = n/2, since weaker shocks are 
the most abundant according to Eq. (|3.9() . Equation l|3.8f) then gives .M* ~ 1.2. 
The maximum depth of the layer is estimated as d ~ 600 pc. With the parameter 
values used to normalize Eq. I|3.12|l and f c \ — 0.25, a = 4.5, we obtain a mean 
separation of canals of L ~ 40 pc, or C ~ 15'. The latter value is reasonably close 
to a by-eye estimate of the average separation of canals in Fig. 3 of Haverkorn et al. 
( 2003 ) of 45'. Thus, the canals observed by Haverkorn et al. H2003|l are compatible 
with a system of shocks producing discontinuities in a Faraday screen. The canals 
in these observations appear to be straighter (or less twisting) than those observed 
at shorter wavelengths by, e.g., Uyamker et al. (1998b) and Gaenslerei al. I|2001[l . 
This might be expected if the straighter canals arise from shocks in the diffuse 
ISM and the twisting canals have a different origin such as differential Faraday 
rotation. The Mach number of the shocks required to produce \AF\ — ir/2 at 
metre wavelengths is quite low, so that they can be inefficient in accelerating 
cosmic rays. Then the expected increase in total synchrotron emissivity at a 
shock is just Ml ~ 1.5. Such features can be difficult to detect in total intensity 
radio maps, especially if their area covering factor is high. 



4 Prospects 

At first sight, the chaotic appearance of radio polarization maps of the diffuse 
Milky Way emission and the difficulty of connecting standard statistical tools, 
such as power spectra, to ISM physics, makes interpretation of the data a daunt- 
ing task. Moreover, changes in the polarization pattern with wavelength depend 
on random quantities, such as the thermal electron density and magnetic field 
fluctuations. Analysis of the ubiquitous depolarization canals provides a practical 
starting point from which we can develop the methods required to make full use of 
the polarization observations. As we have shown here, the origin of canals can be 
determined from observed quantities — by examining how Q and U vary across a 
canal — and their statistical measures are related to the parameters of interstellar 
turbulence. The extension of methods used to study canals to the wider analysis 
of polarized radio emission will lead to better techniques for separating Milky Way 
foregrounds from polarized cosmological signals. 
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